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Abstract 

In Ni-base superalloys with low Al and/or Ti content the precipitation and subsequent coarsening of 
7' particles at intermediate temperatures contribute to the degradation of the mechanical properties 
of the alloy. In the present paper the coarsening process is modelled and the change of the critical 
resolved shear stress of the alloy due to coarsening of the 7' particles is calculated by means 
of statistical analysis of the depinning of a single gliding edge dislocation. It is found that the 
contribution of 7' hardening to the critical resolved shear stress at 973 K reduces to more than 
half of its original value in less than one year. 
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1. Introduction 

Based on their excellent high temperature strength and exceptional resistance to oxidation in 
a wide range of corrosive media, Ni-base superalloys are promising structural materials for next 
generation nuclear power plants. For example, a solid solution strengthened Ni-base superalloy 
(IN617) has been selected as the prime candidate for heat exchanger applications in Generation 
IV power plants [1, 2J. These materials usually contain, among others, a few wt% Al and/or Ti 
solutes. Upon heat treatment precipitation of a fine distribution of 7' particles (Ni 3 Al or Ni 3 Ti) 
is observed with a coherent LI2 crystal structure. These 7' precipitates obstruct the motion of 
dislocations, resulting in high strength. 

After exposure of these alloys to high temperatures for long times, however, the exceptional 
mechanical properties start to degrade. An important reason for this is the gradual coarsening of 
the 7' particles. Since the volume fraction of the 7' phase is nearly constant, the small particles 
disappear and the mean distance between the precipitates monotonously increases with exposure 
time. This implies that the critical resolved shear stress (CRSS) needed for a dislocation to 
overcome the obstacles decreases, resulting in the decrease of the material's strength. In the 
particular case of alloy IN617 (Table [TJ, which is used as a reference throughout the paper, in 
the temperature range between 873 to 973 K, the precipitation of spherical 7' particles with mean 
diameter around 20 nm was observed for a volume fraction of / « 2 — 4% [3[ 0]. After a heat 
treatment at 973 K for 65 600 hours the mean diameter increased to about 200 nm and in parallel, 
a drop in the microhardness was detected [4]. 

The yield stress of a material is usually assumed to be the combination of different strengthening 
mechanisms (such as solute strengthening, work hardening, grain boundary strengthening, etc.). 
In this paper the focus is exclusively on the change in the flow stress due to coarsening of 7' 
precipitates. 
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Element 


Ni 


Cr 


Co 


Mo 


Al 


Fe 


Mn 


Ti 


Si 


C 


Cu 


B 


s 


Fraction (wt%) 


52.2 


22.0 


12.5 


9.0 


1.2 


1.0 


1.0 


0.4 


0.1 


0.1 


0.5 


0.006 


0.015 



Table 1: Nominal chemical composition (wt%) of the Inconel 617 alloy [I]. 



The paper is organized as follows: The 7' coarsening modelled as a diffusion driven Ostwald 
ripening process is described in Section [2j A method based on volume diffusion through the matrix 
is applied. The discrete dislocation dynamics (DDD) model and the procedure to determine the 
CRSS of the alloy is presented in Section [3] The simulation results are then summarized and the 
paper ends with concluding remarks. 

2. Simulation of 7' coarsening 

In the applied model, the precipitates are considered to be hard spheres in a cubic box with 
periodic boundary conditions. To derive the equations of motion governing the coarsening the 
following assumptions are made: (i) the growth is isotropic, (ii) only volume diffusion is present, (iii) 
the diffusion is slow, hence the coarsening is considered quasistatic and (iv) the volume fraction of 
the precipitates is low (maximum few %). The first assumption is usually not valid in 7— 7' systems, 
since the precipitates tend to align in (100) directions, which leads to cubic or other anisotropic 
structures. At low volume fractions and small particles of 7', however, spherical precipitates are 
observed experimentally. In the case of IN617, the 7' particles are spherical even at the latest 
stages of annealing where the mean diameter is around 100-200 nm ECU]. This suggests, that the 
coarsening can be considered isotropic. The differential equations to be solved for the radius of 
the ith precipitate Ri can be then written as [6] 

where n is the total number of particles and the growth rates Bi are calculated from the linear 
equations 

A B = U, (2) 

where A and U are n x n and n x 1 matrices, respectively, and B = (£?i, . . . , B n ). The A matrix 
and the U vector can be constructed at every time step from the coordinates and radii of the 
precipitates - for their actual form (taken from |6J) see the appendix. We note that the variables 
in the above equations are adimensional [5j. 

There are two technical issues that need to be mentioned here. First, the simulations start with 
an initial precipitate number of several ten thousands. The complexity of the solution of Eq. Q 
is 0(n 3 ), and has to be solved at every time step, which is computationally too expensive. The 
complexity can be reduced, however, because there is a screening length A « (R)/f 1/2 for diffusion, 
above which, in principle, no interactions take place between the precipitates [5] (where (R) stands 
for the mean precipitate radius). As such, the simulation cube can be split up into distinct regions, 
and several equations can be solved with lower dimensions simultaneously. The applied method is 
similar to that used by Beenakker [5]. 

The second technical issue is related to the numerical solution of Eq. . For the integration 
we use the 4th order Runge-Kutta scheme with adaptive step size. The latter is necessary, since 
there are slow and fast varying regions during the solution. To demonstrate this, consider the 
event of the disappearance of a small precipitate with index z, i.e. Ri — >• 0. In this case in the ith 
equation of Eq. according to ( |A1[ ) the diagonal element becomes the dominating term in A, 
and the ith equation tends to (l/Ri)Bi = 1/Ri. Thus, in this case Bi — » 1 and the derivative of 
Ri in Eq. diverges. To sum up, at every disappearance of a particle, very small time steps are 
required, while between such events longer time steps are acceptable. 

In practice, the adaptive time step size still does not provide optimal performance. Due to 
the diverging derivatives described above, extremely small timesteps are chosen, but these small 
particles usually have negligible volume (V* ex Rf) and thus they hardly affect the coarsening 
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kinetics of the other particles. To solve this problem, Eq. ([2| was slightly modified by taking U' 
on the right hand side as 

for every z, where for the arbitrary / function the following properties hold: 

1. f(x) > if x > 0, 

2. f{x) —> exponentially as x —> 0, 

3. f(x) — > 1 exponentially as x/R m i n ^> 1, with i? m i n being a characteristic minimum radius 
for the precipitates, 

and (i?)/ := (l/ n ) Si f(Ri)- I n the present case /(#) := [1 — exp(— x/R m i n )] 2 was chosen. It is 
easy to see that if all the radii Ri ^> i? m in, then there is no detectable difference in the equations, 
and when for an i Ri — >• 0, then Bi — >• exponentially, and the derivatives in Eq. ([I]) no longer 
diverge. In addition, the requirement of the constancy of the total volume of the precipitates, 
which is expressed by ^ Bi = [5j, is found not to be violated. To test the validity of the 
outlined approximation simulations were run with different R m i n values as well as without this 
modification, and it was found that (at low enough R m i n values) the introduced error is negligibly 
small. In these calculations, the lattice constant has been taken as the value of R m i n . 

The simulations were performed with a volume fraction of / = 4% and with three different 
cube sizes: L = 1.4 /im, L = 1.75 /im and L = 2.25 /im. Initially, the particles were placed at 
random positions with a lognormal radius distribution, with initial mean radius of Rq — 10 nm 
and relative standard deviation of 20%. Figure [I] shows an example sequence for L — 1.75 /im. 
Note that the initial number of particles was around 20 000 for L = 1.4 /im 50 000 for L = 1.75 
/im and 100 000 for L = 2.25 /im. 
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Figure 1: Coarsening of the 7' particles in the L = 1.75 /im edge box at / = 4% volume content. The mean radius 
(R) increases from 10 nm on the left image to 33 nm on the right. 



The coarsening process has been the subject of extensive theoretical studies. For low (< 10%) 
volume fractions the Lifsitz-Slyozov- Wagner (LSW) theory provides an exact solution for the 
asymptotic (long time) behaviour of the system |8|. In particular, the mean particle radius 
(R) increases as 

(R)(tf-Rl = kt (4) 

at large enough times, where Rq is the mean radius at t = and k is a constant depending 
on the diffusion coefficient, the volume fraction and possibly on the initial spatial distribution of 
the particles. Since there is no available data for diffusion constants in the studied temperature 
range and chemical composition, estimating the coarsening coefficient k relies on derivation from 
experimental observations. Namely, determining (R) is experimentally possible (for example using 
TEM images), from which & exp can be obtained by fitting Eq. Q on such experimental data. This 
was performed on several Ni-base superalloys in the past, mainly with an Al/Ti content in the 
range of 6 — 10 wt.% [9j [TU1 HU H2], and all authors confirmed the coarsening kinetics of Eq. Q. 
More recently, a material with lower 7' volume fraction was investigated by Zhao and co-workers 
|13| . In that case the alloy contained among others, 0.75 wt.% Al and 1.58 wt.% Ti and a similar 
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coarsening behaviour was observed as in the case of IN617 |4J. The measured fc exp followed 



fcexp(r) = k exp 



\RT 



(5) 



with an activation energy of Q = 247 kJ mol -1 and ko = 1.4 x 10 2T nm 3 -h -1 |13] . The coarsening 
rate coefficient k has not been measured for IN617 (to the authors' knowledge), and as such it is 
assumed that it does not differ too much from the value given by Eq. 

It was already mentioned above, that in the simulations dimensionless time units are used. 
However, in order to give a life time assessment for the considered materials, it is necessary to 
link the coarsening process to a real time scale. This can be done using Eq. Q again, which is, 
according to Fig. |2|a), valid in the simulations too. The adimensional simulation time t' can be 
transformed into real time t with a linear transformation t = (/c s im / k exp )t' , where k s [ m is the fitted 
growth rate in simulations, and /c exp is the one measured in the experiments. As such, in 
and[2j the simulation time units are displayed as (k exp /k s i m )t. The linear fit of Eq. Q in Fig 
yields 

k sim = 3.3 x 10~ 2 nm 3 . (6) 

Figure [2^b) shows the evolution of the mean radius (R)(t) itself with the same time units. It can 
clearly be seen that the finite size effects are negligible at these simulation box sizes. 
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Figure 2: (a) The time dependence of the cube mean radius (R) 3 at different simulation parameters. The fitted 
solid line corresponds to Eq. (b) The evolution of the mean radius (R)(t). 



3. Results 

3.1. Simulation set up 

To determine the CRSS of the alloy the threshold shear stress needed for a single edge dislocation 
to continuously penetrate through the 7' precipitates was calculated. Such investigations have 
already been performed several authors on similar materials to the one considered here [T^ [15j UHl 
[13 HH [191 [20] . The main difference is that in the present study, instead of adopting a certain size 
and spatial distribution for the particle arrangement, the data has been obtained from the direct 
simulation of the Ostwald ripening process, described in the previous section. In this work, a 
statistical method [18] is used to obtain the CRSS. Its main features can be summarized as follows. 
The CRSS measured in a single simulation is a probabilistic variable. To obtain the CRSS, which is 
well defined only for systems of infinite size, many smaller systems with different realizations of the 
random precipitate distribution are considered instead of studying a system as large as possible. 
To determine the critical resolved shear stress of an infinite system, the size-dependent transition 
probability Pl is defined as the fraction of realizations where, for a given shear stress and system 
size L the dislocation can traverse the simulation box of size L. Transition probabilities Pl for 
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different system sizes L are well approximated by 



Pl(t) 



1 + erf 



(t - t c )L 
V2a 



(7) 



where L is the system size, a is a fitting parameter for the width of the distribution and r c is 
the critical stress. The form of the distribution function implies, that in the limit of L — >• oo the 
function tends to a step- function at r c , therefore the fitted r c can be considered to be the bulk 
CRSS [18J. 

In the present case, slices were cut from the precipitate configurations of Fig.[I]with appropriate 
plane distances corresponding to the (111) glide planes of the dislocations. The planes were parallel 
to one side of the simulation cube. Then DDD simulations of a single edge dislocation moving in its 
glide plane were performed, as discussed in [17] (for the fundamentals of the DDD method see |21j). 
For the matrix we used fee crystal structure and the following material parameters (corresponding 
to room temperature IN617 |22|): lattice constant of a = 3.66 A, shear modulus of fi = 81 GPa and 
Poisson's ratio of v = 0.3. Two example simulation sequences are shown in Fig. [3] corresponding 
to two different times. 



(a) 



(b) 





Figure 3: DDD simulation of an edge dislocation moving through the particles at constant applied stress and 
L = 1.75 /im. Only the 7' particles intersecting with the glide plane are shown. Sequence (a): t = 10 5 k exp /k s i m , 
sequence (b): t = 10 6 fc e xp/^ s im- 



It has to be noted, that while the coarsening process is simulated in the temperature range 
of 873—1073 K, the CRSS determination is performed at room temperature by averaging over 90 
realizations. This means that only dislocation glide is assumed, and temperature induced processes, 
such as dislocation climb or cross slip, are neglected. 

In Ni-base superalloys, the transition temperature between the Orowan process and other high 
temperature bypassing mechanisms is around 1073 K [23]. Below this value the relative drop in 
the yield stress can be attributed to the same relative drop in the elastic constants, which suggests 
that the bypassing mechanism is not changed. Our DDD model is, therefore, valid also in the 
temperature range of 873—1073 K, only the elastic constants should be modified. To extend these 
studies to even higher temperatures the inclusion of dislocation climb and cross slip into the model 
is the subject of further work. 

3.2. Determination of the CRSS 

The measured cumulative distribution functions of the CRSS of individual 2D slices for the 
L = 1.75 /im sample are plotted in Fig. [4] at different times. Each of the curves is constructed from 
~ 90 simulations. 
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Figure 4: The cumulative distribution function of the critical shear stress values. The solid and dashed colour lines 
correspond to a system size of L = 1.75 fim and L = 2.25 /im, respectively. The thin black lines are the fitted 
functions P L of Eq. Q for L = 1.75 



The determined CRSS as a function of time is shown in Figure |5j In the measured range the 
CRSS drops more than 50%, while the mean radius of the particles increases from 10 nm to about 
35 nm. To relate the dimensionless timescale to the real one, Eq. ([5| is implemented with T = 973 
K. In this case /c ex p — 4.22 nm -h , and k s i m 

Aexp = 7.82 x 10 -3 h. The top scale of Fig. [5] 
corresponds to this transformation, and shows that the drop of more than 50% in the CRSS due to 
coarsening of precipitates happens in less than one year, even at this moderate temperature. For 
a more precise estimation for IN617, the experimental measurements of the coefficient /c exp would 
be necessary. 
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Figure 5: The evolution of the CRSS. The upper timescale assumes T = 973 K. For comparison, the (R)(t) function 
was plotted again, as well as the theoretical expression of Eq. (pH) (green line) for the CRSS. 



Several analytical expressions exist for the CRSS due to a distribution of precipitates bypassed 
by the Orowan mechanism [24] . These are based on assumptions on the spatial and radial distribu- 
tion of the particles. In our case, the coarsening starts from spatially random initial configurations 
with a lognormal size distribution (see Section [2| . During the Ostwald ripening, both spatial 
correlations of the particles evolve and the size distribution also tends to an asymptotic function, 
described within the LSW theory [7, 8J. To estimate the effect of these phenomena on the CRSS 
our results are compared with the formula derived by Nembach based on the work of Bacon and 
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co-workers (23 [2§Q 

r p ((R)) = ^^H^n(R)/b). (8) 

In the formula co>l and co>h are numeric constants depending on the geometric properties of the 
precipitates. Here the values of co>l = 5.69 and ujh = 1.22 were adopted, which were derived for 
monodisperse particles [24]. In Fig. [5] the solid line corresponds to Eq. (|8|, without any fitting 
parameters. The remarkably good agreement suggests that even if the spatial correlations affect 
the CRSS, the dependence must be weak or only effective in case of particle clustering and particle 
free areas. 



4. Conclusions 

A new method was established to study the effect of j' coarsening on the CRSS of Ni-base 
superalloys with a low volume fraction of 7' phase. The assumptions taken in this model are 
that the coarsening is driven by bulk diffusion through the matrix (Ostwald ripening), and that 
the volume fraction of the precipitates remains constant during the growth. The timescale of 
the growth process was determined by fitting to experimental observations. The critical resolved 
shear stress of the material was measured at different stages of the coarsening with a DDD method 
coupled with a statistical analysis. It was found that the spatial inhomogeneities introduced during 
the coarsening do not have significant effect on the CRSS and that in an IN617 alloy at 973 K the 
CRSS contribution of the 7' phase decreased by more than 50% of its initial value in less than 1 
year. 



Appendix 

The form of the matrices introduced in Eq. ([2| |6J: 



/ 71 r> 71 r> 

/ J l_ R± _i l_ R± 

Ri n{R) ^ rji ri2 n{R) ^ r j2 



1 

T21 



1 yy Rj 



i{R) 



1 V ^± 

% (R) ^ r n 



1 

R2 



1 

r n 2 



1 R j 



j/2 



Rj_ 



1 

Tin 



1 

r 2n 



1 

R n 



_JL_ f *l\ 

1 Rj_ 
^ r 3n 



n{R) 



(R) ^ r 



R± 



where denotes the distance between the ith and jth particles. 
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